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Abstract 

The a-bosonic properties such as single-a orbits and occupation numbers in J'^=0^, 2"*", 1~ 
and 3^ states of ^^C around the 3a threshold are investigated with the semi-microscopic 3a 
cluster model. As in other studies, we found that the 0^ (2^) state has dilute-3a-condensate-like 
structure in which the a particle is occupied in the single S {D) orbit with about 70 % (80 %) 
probability. The radial behaviors of the single-a orbits as well as the occupation numbers are 
discussed in detail in comparison with those for the 0^ and 2f states together with the lj~ and 3j~ 
states. 

PACS numbers: 21.10.Dr, 21.10.Gv, 21.60.Gx, 03.75.Hh 
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I. INTRODUCTION 



Four-nucleon correlations o 
microscopic a-cluster model 



" the a-cluster-type play an important role in nuclei. The 



has succeeded in describing the structure of many 
states in light nuclei, in particular, around the threshold energy of breakup into constituent 
clusters. As for ^^C, detailed analyses were made by several authors with the microscopic 
3a cluster model about 25 years ago. The 3a GCM (generator coordinate method) 6] and 
3a RGM (resonating group method) calculations showed that the second O"*" state of ^^C 
{Ex=7.65 MeV), located at _E3a=0.38 MeV above the 3a threshold, has a loosely coupled 3a 
structure, although the ground state is a shell-model-like compact state. On the other hand, 
special attention has been paid to an a-type condensation in symmetric nuclear matter, 
analogue to the Bose-Einstein condensation for finite number of dilute bosonic atoms such 
as »'Rb o. ^»Na at ve.y low te,„pe.at.e whe.e al. ato^s occup. the lowest . o.b.t fl. 
Several authors 0, 3| showed the possibility of such a-particle condensation in low-density 
nuclear matter, although the ordinary pairing correlation can prevail at higher density. This 
result suggests that dilute a condensate states in finite nuclear system may exist in excited 
states weakly interacting gas of a particles. 

Recently, a new a-cluster wave function was proposed which is of the Na particle con- 



densate type 



\<^Na) = {C^r\y^c), (1) 

(r-i ■ ■ ■ r;v|$iv.) oc A |e-K^?+- +n)0(tt,) . . . , (2) 

where is the a-particle creation operator, 0(a) denotes the internal wave function of 
the a cluster, is the center-of-mass coordinate of the i-th a cluster, and A presents 
the antisymmetrizer among the nucleons belonging to different a clusters. The important 



uster is of S- 
indicated 



characteristic of the wave function is that the center-of-mass motion of each a c 
wave type. Applications of the condensate- type wave function to ^^C and ^^O 
that the second 0+ state of (E^=7.65 MeV) and fifth 0+ state of ^^O (^^=14.0 MeV), 
around the 3a and 4a threshold, respectively, are conjectured to be dilute A^a condensate 
states, which are quite similar to the Bose-Einstein condensation of bosonic atoms at very 
low temperature. The calculated nuclear radii for both of those states are about 4 fm, 
significantly larger than that for the ground state (about 2.5 fm). As for ^^C, a detailed 
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analysis with a deformed alpha condensate wave function, slightly different from the spherical 
one in Eq. (0), was performed to investigate the structure of the 0^ and 0^ states. It was 
found that each of the 0^ wave functions obtained by the 3a GCM and RGM calculations 
has a large squared overlap value of more than 90% with the single 3a; condensate wave 
function. 

The above-mentioned results for ^^C and ^^O lead us to the further intriguing issue that 
dilute a-cluster states with J'^ = 0^ near the Na threshold may exist in other heavier 
AN self-conjugate nuclei. The Gross-Pitaevskii-equation approach [l^ is useful to explore 
such dilute multi a systems, because this equation [l^l, based on mean field theory, has 
succeeded in describing the structure of the Bose-Einstein condensation for dilute neutral 
atomic systems, for example, ^"^Rb or ^^Na, at very low temperature, trapped by an external 
magnetic field. Q The present authors [l^ applied the Gross-Pitaevskii equation to self- 
conjugate AN nuclei. They found that 1) there exists a critical number of a bosons that the 
dilute Na system can sustain as a self-bound nucleus, and 2) the Coulomb-potential barrier 
plays an important role to confine such dilute A^a-particle condensate states. 

It is interesting to explore also the possibility of the a condensate states with non-zero 
angular momentum in ^^C. The old theoretical calculations based on the microscopic 3a 
cluster model {4,0,01 suggested the existence of a 22 state of ^^C at around E^a ~ 3 MeV 
above the 3a threshold, the structure of which is similar to the O2" state except for the 
angular momentum. Quite recently the 22 state was observed at E^a = 2.6 ± 0.3 MeV with 
the alpha decay width = 1.0 ± 0.3 MeV. [3| The a condensate- type wave function 



with axially deformation jlq] was applied to study the structure of the 2 J st ate with help 
of the method of ACCC (analytic continuation in the coupling constant) jl7|. They found 
that the 22 state has a large overlap with the single condensate wave function of 3a gas-like 
structure, the squared value of which amounts to about 88%. This result implies that the 
22 state has a similar structure as the 02" state, namely, dilute 3a condensation. 

Here, it is an intriguing problem to discuss whether the 0^ and 2^ states of ^^C are 
ideal dilute a condensates or not. The condensate-type a-cluster wave function in Eq. Q 
has succeeded in describing the 0^ state of ^^C. This result, however, does not necessarily 
mean that the 0^ state of ^^C is an ideal a-condensate state. If the 02" state of ^^C is 
an ideal dilute a-condensate, the single a-particle orbit in the state should be of the zero- 
node long-ranged S-wave type with an occupation probability of 100 %, as suggested from 
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the Gross-Pitaevskii-equation approach The antisymmetrizer A in Eq. Q generally 

perturbs the single a motion in the nucleus, and one should remember that the condensate- 
type wave function can also describe the shell-model-like compact structure of the ground 
state of ^^C. The effect of the antisymmetrizer should have a close relation to the rms radius 
of the nucleus or the distance between 2a clusters in a nucleus. The ideal 3a condensate state 
is expected to be realized if the distance between two arbitrary a clusters is large enough 
so that the effect of the antisymmetrizer can be neglected. The calculated nuclear radius 
for the 0^ state of ^^C, about 4 fm h|, suggests that the action of the antisymmetrizer 
is weakened significantly in that state. In order to give more decisive theoretical evidence 
that the 0^ state of ^^C as well as the 2^ state has dilute 3a condensation structure, it is 
needed to study quantitatively the bosonic properties such as single a-particle orbits and 
corresponding occupation probabilities, starting from the microscopic wave function. 

The first attempt to derive the a-boson properties for 0''" states in ^^C from a microscopic 
model was performed in Ref. where the 3a RGM equation was solved in terms of the 
correlated Gaussian basis with the stochastic variational method. Although they formulated 
a derivation of the 3a boson wave function starting from the microscopic 3a wave function, 
the a bosonic properties of ^^C were studied not with the 3a bosonic wave function but 
with the normalized spectroscopic amplitude, because the derivation of the 3a boson wave 
function is numerically difficult due to the non-local properties of the norm kernel. Although 
the normalized spectroscopic amplitude seems to be a good approximation for the boson wave 
function in the region where the effect of the antisymmetrizer is negligible, the approximation 
becomes worse when the spatial overlap of the 3a clusters becomes larger. It is requested 
to demonstrate quantitatively how good the approximation is for the 0^ state within their 
framework. 

The purposes in the present paper are twofold. First we study the bosonic properties 
such as single a-particle orbits and occupation probabilities for the O"*" and 2"*" states in 
^^C with direct use of the wave function obtained by the 3a OCM (orthogonality condition 
model) The OCM is a semi-microscopic model and a simple version of the RGM, taking 
into account properly the antisymmetrization among nucleons, which successfully describes 



the structure of ^^C. 
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I The second purpose is to explore the possibility of the 
dilute 3a condensation with negative parity within the present framework. The 3j~ 
state of at E^a = 2.37 (3.57) MeV above the 3a threshold appears at the same energy 



4 



region as that for the 0^ (£"^=0.38 MeV) and 2^ (2.6 MeV) states. According to the old 
theoretical study based on the 3a GCM and RGM calculations 0, Q|, the nuclear radius 
of the 3~ state is intermediate between a compact shell-model-like state {Of) and a loosely 
coupled 3a cluster state (0^), while the 1~ state has a radius only a little smaller than that 
of the 0^ state. Thus, it is quite interesting to study the bosonic properties for the negative 
parity states, as well. 

In order to clarify the relationship between the Na boson wave function and the one of 
Na OCM, we first outline a way of mapping of the microscopic Na wave function onto the 
Na boson wave function, and derive the equation of motion for Na bosons from the Na 
RGM equation in Sec. II. The Na OCM equation is illustrated as an approximation of the 
Na boson equation. The Na OCM wave function, thus, has bosonic properties. The 3a 
OCM equation is solved properly with modern numerical techniques. The calculated single 
a-particle orbits and occupation probabilities in ^^C are discussed in Sec. III. Finally, we 
give the summary in Sec. IV. 
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II. FORMULATION 



A way of mapping of the microscopic Na wave function onto the Na boson wave function 
is illustrated in order to derive the equation of motion for the Na bosons. The Na OCM 
equation is given as an approximation of the equation of motion for Na bosons. We formulate 
the evaluation of the single-a orbits and occupation numbers from the Na OCM wave 
function together with other physical quantities. Finally, we give an outline of how to solve 
the 3q; OCM equation for ^^C with a phenomenological aa potential. 



A. Mapping of the fermionic Na wave function onto a Na boson wave function 

In the microscopic Na cluster model, the total wave function, \E'j , with the total angular 
momentum J is given as 



¥P =A{1[ <p'^fxj (r) \ , (3) 



N 



da¥f\a)x{a), (4) 

where (j)^^*^ denotes the intrinsic wave function of the a particle with the simple (Os)^ 
shell-model configuration, and xj represents the relative wave function with a set of the 
relative (Jacobi) coordinates, r = {ri, 7-2, ■ ■ ■ , Ttv-i}, with respect to the cm. of a clusters. 
The antisymmetrization among 4N nucleons is properly taken into account in terms of the 
operator A. The function '^^{a) in Eq. (jlj) is defined as 

{N N-i ^ 
i=i j=i J 

which describes the a-cluster state located at the relative positions specified by the Jacobi 
parameter coordinate a = {ai, a2, ■ ■ ■ , a^^i}- 

The total Hamiltonian for 4N fermions is given as 

AN AN 
1=1 i<j=l 

where ti and Tcm denote, respectively, the kinetic energy operator of the i-th nucleon and of 
the center-of-mass of the total system. The nucleon-nucleon interaction (Coulomb interac- 
tion) between the i-th and j-th nucleons is expressed as Vij 
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The Schrodinger equation for the fermionic Na system is 

Hmf^ = E¥p. (7) 

Substituting the total wave function of Eq. (jD) into Eq. ((Tj), we obtain the equation of motion 
for the relative wave function xji 

Jda'{H{a,a')-EN{a,a')}xj{a') = 0, or {n-Ar)xj = 0, (8) 

which is called the RGM (resonating group method) equation |5|]. The Hamiltonian and 
norm kernels, H{a,a') and N{a,a'), are defined as 

= {<\a) I I ^ I I (a')). (9) 

Recalling the normalization condition {^j \ ) = 1 for the total wave function in 
Eq. Q, the normalization of xj in Eq. (jH} is given by 

J dada'x*Aa')N{a, a')xAa') = 1. (10) 

( Fi) 

This suggests that an Na boson wave function <l>j corresponding to the fermionic wave 
function in Eq. Q should be taken to be 

<!>f{a)=Af'/\j = J da'N'/\a,a')xAa'), J da\<i>f\a)\' = 1, (11) 

where iV^/^ is defined as 

I da"N'/\a, a")N^/^{a", a') = N{a, a'). (12) 

It is noted that the boson wave function ^[f \a) has only the Jacobi coordinates a = 
{ai, a2, ■ ■ ■ , <2Af-i} of the Na system and the internal coordinates in the a cluster are in- 
tegrated out completely. In addition, $j (a) is totally symmetric for any two-a-cluster 
exchange. Thus, it has bosonic property. From the RGM equation (jHJ, should satisfy 
the following equation 

(Ar-V27^A/--V2 _ $(f ) = 0, (13) 

Thus, we can interpret J\f~^^'^'HN'^^^'^ as the Na boson Hamiltonian, and Equation ()13|) is 
the equation of motion for the Na boson wave function. 



If the eigenvalue problem for the norm kernel 

Mux = J da'N{a, a')ux{a') = Xux{a) (14) 
is solved, the boson wave function ^j^\a) in Eq. (fTT|) . then, is obtained as 

^^f\a) = Y.^Ma){u,\xj). (15) 

A 

The eigenvalue of the norm kernel. A, is non-negative, and the eigenfunction u\ with A = 
is called the Pauli-forbidden state, which satisfies the condition ^{Hili 0af*^'WA} = 0. The 

( 

boson wave function $j , thus, has no component of the Pauli-forbidden state, 

{ux I $f ^) = for Ma with A = 0. (16) 



B. Na orthogonality condition model (OCM) with bosonic properties 

In the previous section, we mapped the fermionic wave function onto the Na boson wave 
function ^^f^ within the framework of the resonating group method (RGM). The boson 
wave function has the following properties: 1) $j is totally symmetric for any 2a-particle 
exchange, 2) satisfies the equation motion in Eq. (fT^ . and 3) is orthogonal to the 
Pauli forbidden state (see Eq. ()16p). In order to obtain the boson wave function, we need to 
solve the RGM equation (jH)) and the eigenvalue equation of the norm kernel (jl4p or to solve 
directly the equation of motion (fT!?|l . Solving the eigenvalue equation of the norm kernel, 
however, is difficult in general even for the 3a case. In addition, computational problems 
are encountered for solving the Na RGM equation (jSI) for > 4. Thus, it is requested to 
use more feasible frameworks for the study of the bosonic properties and the amount of a 
condensation for the Na system. In the present study, we take the orthogonality condition 
model (OCM) as one of the more feasible frameworks. The OCM scheme, which is an 
approximation to the RGM, is known to describe nicely the structure of low-lying states 



in light nuclei |^, Q, IZll |22 1 . The essential properties of the Na boson wave function 
$j , as mentioned above, can be taken into account in OCM in a simple manner. We will 
demonstrate this briefly. 

In OCM, the a cluster is treated as a point-like particle. We approximate the Na boson 
Hamiltonian (non-local) in Eq. (jl3|) by an effective (local) one H^'~"~^^\ 

x-'/^nx-'/^ ~ (17) 



N N N 

^(OCM) ^Y.T^-nn.+ E V.^^^ (t, j) + ^S^l^, J, k) , (18) 

i=l i<j=l i<j<k=l 

where Tj denotes the kinetic energy of the «-th a cluster, and the center-of-mass kinetic 
energy Tcm is subtracted from the Hamiltonian. The effective 2a and 3a; potentials are pre- 
sented as V2a^ and V^^-^ , respectively. Referring to the RGM framework in Eqs. (fT^ ^(fTH |) . 
the equation of the relative motions for the Na particles with H^'~"^^\ called the OCM 
equation, is expressed as 

l^(OCM) _e}^j = 0, (19) 
{up I $j) = 0, (20) 
where up denotes the Pauli- forbidden state, satisfying the following condition 

N 

Am<l>^fuF} = 0. (21) 

i=l 

The orthogonality condition in Eq. (PUI) corresponds to Eq. (fT^ in Sec. IIA. The bosonic 
property of the wave function $j can be taken into account by symmetrizing the wave 
function with respect to any 2a-particle exchange, 

$j = 5$j(l,2,---,iV), (22) 

where S denotes the symmetrization operator, S = {I / VN\) J2v ^ where the sum runs 
over all permutations V for the Na particles. It is noted that the complete overlapped state 
of the 3a particles is forbidden within the present framework because of the Pauli-blocking 
effect in Eq. (j2DI), although we take into account the bosonic statistics for the constituent a 
particles. In the next section, we will demonstrate i) how to solve the OCM equation and 
ii) what kind of effective a-a potential we should choose in for the 3a case of ^^C. 

Here, it is useful to define various quantities characterizing the structure of the Na 
system with use of the Na boson wave function $j obtained by solving the OCM equation 
in Eqs. IH^ and (pUjl . The single a-particle density is defined as 

N 

pir) = {<i>j\Y.^ir-r\-^)\<^j), (23) 

i=l 

where r •'^"^^ is the coordinate of the i-th a particle measured from the center-of-mass coor- 
dinate of the total system. The nuclear root-mean-square radius is given as 



{r%) = ^{rl) + l.7l\ (24) 
W)= fdrr'pir), (25) 



where we take into account the finite size effect of the a particle. The correlation functions 
with respect to the a-a relative coordinate Vaa as well as the relative coordinate between 
one of the a particles and the remaining [N — V)a system roL-{N-\)a are given as 

/aa(r) = ($j|5(r-r,J|<l>j), (26) 
/a-(iv-i)a(r) = ($,/ I 8{r - | (27) 

where the way of choosing the coordinates, Vaa and ra-{N-\)a-, is arbitrary in the set of 
Jacobi coordinates of the Na particles because of the totally symmetrization for $j. The 
root-mean-square (rms) distances with respect to r^a and ra_-{N-\)cx are, respectively, given 
by 

1/2 



(rL)= k^AKA^j) , (28) 



1/2 



(rL(^_i)J = I I . (29) 

The reduced width amplitude for the a-{N — l)a part is defined as 

yiLjir) = r X ([r^(r)0, ((iV - l)a)]j \ $,) , (30) 

where r denotes the relative coordinate between the a particle and the (A^ — l)a nucleus, 
and (pe {{N — l)a) represents the wave function of the (A^ — l)a nucleus with total angular 
momentum i which is obtained by solving the OCM equation for the (A^ — l)a system. The 
integration in Eq. is done over all of the relative (Jacobi) coordinates for the Na system 
except for the radial part of r. 

In order to discuss the bosonic properties such as the degree of a condensation in a 
nucleus, it is needed to extract information on the single a-particle orbits and its occupation 
probabilities in the nucleus from the total wave function $j. The one-particle density matrix 
for the Na system is very useful for this Defining the one-particle density operator as 

N 

^ir,r') = Y: I 5ir\^-^ - r')) (5(rJ'^"^) - r) |, (31) 

2 = 1 

then, the single a-particle density matrix is given as 

p(r,r') = {<l>j\V{r,r')\<^j), (32) 
= Nx{^j I 5(ri^"^) - r'))(5(ri^"^) - r) | (33) 
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where is the same as that in Eq. ()23|) . It is noted that the diagonal matrix element 

reduces to the single a-particle density in Eq. (P^: p(r, r' = r) = p{r). The single a-particle 
orbit and its occupation number in the nucleus can be evaluated by solving the eigenvalue 
equation of the one-particle density matrix 

j dr'p{r,r')ip^ir') = p^^{r), (34) 

where the eigenvalue /i presents the occupation number. The eigenfunction denotes the 
single-a orbital wave function in a nucleus with the argument of the intrinsic coordinate 
(^^(cm)-j q£ arbitrary a particle in a nucleus measured from the center-of-mass coordinate. 
The ratio fi/N represents the occupation probability of an a particle in the orbit (f^^. The 
spectrum of the occupation probabilities offers important information about the occupancy 
of the single a-particle orbit in a nucleus. If each of the Na particles in an A^a-boson state 
is occupied by only one orbit, the occupation probability for this orbit becomes 100%. 

The *Be (2a) system is a good example to demonstrate the characteristic of the single-a 
orbital wave function. From the definition of Eqs. and ()34|1 . the L^-wave single-a orbit 
in the ®Be( J'^) state corresponds to the relative wave function (which is obtained by solving 
the 2a OCM equation with J {= La) in Eqs. f|TT?|l and (f^ ). scaling to 1/2 with respect to 
the relative coordinate between the 2a clusters. Then, the occupation probability becomes 
exactly (mathematically) 100% for any L^j-value. 

The radial behavior of the L^-wave single-a orbit, ffj.{r'^"^^), in Eq. ()34j) generally has a 
close relationship with that of the reduced width amplitude, y£Lj{ra-{N-i)a), in Eq. (jHS]). 
This is due to the fact that both represent the behavior of the single-a-particle motion in a 
nucleus in which all degrees of freedom of the other a particles are integrated out, and r'^"^^ 
is given as r^f"") = ^ x ra-{N-i)a- 

The momentum distribution of the single a particle is also helpful for the study of a 
condensation in a nucleus |lSj. It is defined as a double Fourier transformation of the 
one-particle density matrix 

p{k) = J dr'drj^^pir, r')j^,, (35) 

J dkp{k) = 1, (36) 

It is reminded that p{k) would have a ^-function like peak around k = for an ideal dilute 
condensed state of infinite size. 
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C. 3a OCM for ^^C 



In the previous section, we outlined the Na orthogonahty condition model (OCM) and 
discussed how to extract the properties and the amount of a condensation in the Na system. 
Here, we will apply the OCM framework to the 3q; system of ^^C. 

The total wave function of ^^C (the total angular momentum J) within the frame of the 
3q; OCM is presented as 

where $j ' denotes the relative wave function of the 3a system with the Jacobi-coordinate 
system shown in Fig. Hl^a), and others are_ self-explanatorv. In the present study, $ j(^^C) is 
expanded in terms of the Gaussian basis 



<^'irc) = EE^c(^,/i)$r(^,/x), (38) 
= $f + + $f (39) 

<|.p)(z.,^) = [cp,ir,,,u)Mrk,f^)]j, (40) 
ipe{r,i^) = iV,(z/)/exp(-z/r2)F,(f), (41) 

where Ni is the normalization factor, and r^j (r^) denotes the relative coordinate between 
the i- and j^th a particle (the k-th a particle and the center-of-mass coordinate of the i-th 
and j-th a particle). The angular momentum channel is presented as c = {i,L)j, where i 
(L) denotes the relative orbital angular momentum between 2a clusters (the center-of-mass 
for the 2a clusters and the third a). The Gaussian parameter u is taken to be of geometrical 
progression, 

z/„ = 1/bl, hn = &mina"~\ n = 1 ~ nmax- (42) 

It is noted that the prescription is found to be very useful in optimizing the ranges with a 

small number of free parameters with high accuracy [2^ . 

The total Hamiltonian for the 3a system is presented as 
3 3 

n = Y.T,~T,m+ E [V2a{nj) + V^^'^'in,)] + V^a{ri2, r23, rsi) + Vpauli, (43) 
1=1 i<j=l 

where T,, and V^a stand for the kinetic energy operator for the i-th a particle, phe- 
nomenological 2a and 3a potentials, respectively, and V^*^°"' is the Coulomb potential be- 
tween 2a particles. The center-of-mass kinetic energy is subtracted from the Hamiltonian. 
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The Pauli-blo eking operator Vpauii is represented as 

^Pauii = lim AOpauii, (44) 

A— >oo 

Opauli = Y1 \Une{rij){UnE{rij)\, (45) 

2n+e<4.,e=even (ij)={12),{23),{31) 

which removes the Pauh forbidden states, Mqo, Mio and U20, between any two a particles from 
the 3a model space. The Gaussian size parameter of the nucleon (Os) orbit in the a cluster 
is taken to be bj^i = 1.358 fm, which reproduces the size of the a particle in free space. In 
the present study, we take the harmonic oscillator wave functions as the Pauli forbidden 
states. The eigenenergy E of ^^C and coefficients in Eq. (jHHj) are obtained in terms of 
the variational principle, 

6 [{<^j \H-E\ <l>j)] = 0. (46) 

We use an effective 2a potential which reproduces the observed aa scattering phase shifts 
{S-, D- and G-waves) and the resonant ground-state energy within the 2a OCM framework. 
The effective 2a potential and Coulomb potential, and Vg^""', are constructed with the 
folding procedure, where we fold the modified Hasegawa-Nagata effective NN interaction 
(MHN) and the pp Coulomb potential with the a-cluster density. Also the strength of the 
second-range triplet-odd part in MHN is modified so as to reproduce the 2a scattering phase 
shifts. 

Only using the effective 2a potential leads to a significant overbinding energy for the 
ground state of ^^C within the frame of the 3a; OCM. Thus, we introduce an effective, 
repulsive, 3a potential, V30, in addition to the 2a potential, 

V,^ = Vo exp [-P {rl, + rj, + J] , (47) 

where Vij denotes the relative coordinate between the i- and j-th a particles, and Vq and (3 
are taken to be Vq = 87.5 MeV and f3 = 0.15 fm~^. Including the 3a potential, the energy 
of the ground state of ^^C is reproduced with respect to the 3a threshold, together with the 
nuclear radius (see Sec. III). 

Single-a orbits and corresponding occupation probabilities for O"*", 2+, 1~, and 3~ states 
of ^^C are investigated by solving the eigenvalue equation of the single a-particle density 
matrix in Eqs. (jHHj) and [see Sec. 11(b)]. They will lead to a deep understanding about 
the structure of ^^C. 
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In the present investigation, we make a further structure stu dy f or the 0^ states of 
^^C, because they have very intriguing features. According to Ref. the Of state has a 
compact shell-model-hke state, while the 0^ one is conjectured to have a dilute 3a condensate 
structure, the nuclear radius of which is 4.3 fm, much larger than that of the ground 0^'" state 
(2.48 fm). Thus, it is interesting to see in detail the structure change of the 0"*" state of ^^C by 
taking the nuclear radius as a parameter. We investigate the dependence of the occupation 
probabilities and radial behaviors of the single a-particle orbits in the 0"^ state on its rms 
radius within the 3a OCM framework. The results will give us helpful understanding about 
the structure of ^^C. The procedure of evaluating them is formulated hereafter. 

First, we consider a Pauli-principle respecting 3a OCM basis wave function. For the 
purpose, the eigenvalue problem for the Pauli operator in Eq. (j45|) is solved to obtain the 
Pauli forbidden state in the 3a OCM model space 

OpauiilGp") = AplGp"), (48) 

where Xp denotes the eigenvalue for the eigenfunction IGj?). The Pauli operator, then, is 
expressed as 

Opauii = El^p)^^(G'p I- (49) 
p 

If Xp is non-zero, its eigenfunction corresponds to the Pauli forbidden state. In the present 
study, the eigenvalue problem is solved with use of the 3a OCM basis in Eq. (j39|) . Then, 
the Pauli-principle respecting OCM basis wave function is given by 



<l>f(z/,/x) = iV,(z/,/i) 



(50) 



where denotes the normalization factor with the angular momentum channel c = {i, L)j, 
and $j?"(z/, /i) is given in Eq. 

According to the results of the 3a OCM calculation (see Sec. IIIA), the ground state (O^") 
and second 0^ states of ^^C have the equilateral triangle configuration of the 3a clusters. In 
addition, it is found that the single- angular-momentum-channel calculation with c = {iL)j = 
(00)o gives a good approximation to the results of the full coupled-channel calculation. 
Thus, only the single angular momentum channel c = {iL)j = (00) o is taken in the present 
calculation, and the equilateral triangle configuration is assumed for the basis wave function. 
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The latter can be realized easily by putting the condition u = fi in Eq. (ISUp. Then, the rms 
radius of the Pauli-principle respecting 3a OCM basis wave function is evaluated as 



1/2 



(51) 



1=1 / 
which depends on z/, and where we take into account the finite size of the a cluster. The 
energy of the 3a system is given by 



E^M = ($Jo)o(^'/^ = ^)\n\<l>t^,)^{u,f, = v)) . (52) 

where 7i denotes the total Hamiltonian of the 3a system in which we subtract the Pauli- 
blocking operator Vpauii from the 3a OCM Hamiltonian Ti in Eq. (jlHj) . The single-a density 
matrix is given as 

p^{ry) = (l-foVl^./^ = ^) I E I ^(^f^^ - r')){6{r\^-^ - r) \ ^l^,^,{v,^^ = v)). (53) 

i=l 

The single a-particle orbit and its occupation number in the basis wave function are obtained 
by solving the eigenvalue equation of the single-a density matrix, 

j dr'p^{r,r')ip.n{r') = r](pr,{r), (54) 

where the eigenfunction 99^ denotes the a-particle orbit with the occupation number 77 (eigen- 
value). Thus, we can study the dependence of the occupancy of the single a-particle orbits 
in the 0"*" state of ^^C on its nuclear radius by choosing the parameter value u so as to 
reproduce a given nuclear radius. 
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III. RESULTS AND DISCUSSION 



A. 0^ and O2 states 

Table HI presents the results for the energy, measured from the 3a threshold, and nuclear 
radii for the ground (Of) and excited states (O2 ) of ^^C. The energy for the ground state 
is reproduced well, and the corresponding nuclear radius, 2.44 fm, is in good agreement 
with the experimental charge radius (2.4829 ± 0.019 fm) with an error of about 2%. On the 



other hand, the rms distance between 2a clusters in the Oj" state is y {r'^)aa=3.02 fm (see 
Table H} , which is larger than that between the center-of-mass of the 2a clusters and the 
third a cluster, ^ (r2)^„2a=2.61 fm. Then, the square of the ratio, ^ {r'^)a-2a/ \J {'^'^)aa , is 
about 3/4. The results mean that the ground state has an equilateral-triangle- like intrinsic 
shape. Figure El shows the density distribution of the a particle for the O^" state of ^^C. We 
see a prominent peak at r ~ 2 fm, which demonstrates clearly the shell-model-like compact 
structure of the ground state of ^^C. 

As for the 02" state, the energy measured from the 3a threshold is -E3q,=0.86 MeV 
{E^=8.13 MeV), which agrees well with the experimental data -£^3^^=0.38 MeV (i?|^'^=7.65 
MeV). The calculated nuclear radius for the 02^ state is as large as 4.31 fm (see Table HJ. 
This means that the state has a dilute 3a structure, although our nuclear radius is a little 
large. that in Ref. Q. The dens.ty d.t.bution of the a particle for the .ate Is 
illustrated in Fig. El In comparison with that for the ground state, we can easily recognize 
the dilute structure of the 02" state, which is in contrast with the compact structure of the 
ground state. 

The difference between the structures of the 0^^ and 0^ states can be also seen in the 
radial behavior of the correlation functions, faa and fa-2a, with respect to the a-a and 
a-2a relative coordinates, respectively, of Eqs. and They are illustrated in Fig. O 
Reflecting the compact structure of the Of state, both of faa and fa-2a have prominent 
peaks at r ~ 2.6 fm and 2.5 fm, respectively, and extend to r ~ 5 fm, while those for the 
02" state show bump structures with peaks at r ~4 fm and r ~5 fm, respectively, and have 
a long tail up to r ~ 15 fm. 

It is instructive to study the single a-particle orbits (eigenfunctions) and occupation 
numbers (eigenvalues) of the one-body density matrix in Eq. ((221) • The results of the di- 



16 



agonalization of the latter are shown in Table |n] together with the occupation probability 
defined as the occupation number divided by the number of a particles. As for the ground 
state, the occupation probabilities spread out over S", D and G waves, but they are concen- 
trated to the first orbits, Si, Dx and G\ orbits, respectively, where denotes the fc-th orbit 
for the L wave. The occupation probabilities are about 30 % for all orbits. This result is ex- 
pected from the fact that the ground-state wave function is of nuclear SU(3)-like character, 
SU(3)[/](A/i)j. = [444](04)o+ with quanta Q=8, where the SU(3) bases with Q < 8 corre- 
spond to the Pauli-forbidden states. Since the SU(3)[444](04)o+ state is the eigenfunction 
of the 3a RGM norm kernel in Eq. ()14j) . it can be regarded as the 3q; boson wave function 
with Q = 8, see Eq. (fT^ . The state is described as 

|[444](04))o+ = an;iVL|(n£)(iVL)), 

nlNL 

where \{nt){NV)) presents the basis function \unt{r 2a)u ^i^{ra-2a)) with 2n + £+2A^+L = 8, 
and Unt (unl) denotes the harmonic oscillator wave function with the number of nodes n 
(N) and orbital angular momentum i (L) referring to the coordinate vector {ra-2a) 
between 2a clusters (between the center-of-mass for the 2a clusters and the third a cluster). 
Let us define Lq, as the orbital angular momentum of a single-a orbit. Then, L in Eq. ()55p 
corresponds to La, because La is defined as the orbital angular momentum with respect to 
^(cm)^ coordinate vector of the a particle measured from the center-of-mass coordinate of 
^^C (see Eq. fj23|) ). which is parallel to ra-2a (t'q'"^ = ^fa-2a)- From the definition of the 
one-body density matrix in Eq. ()33|1 . the single-a orbits and occupation probabilities for the 
SU(3) state in Eq. ^ are given as follows: 64/255 ~ 28% for 5-orbit, 80/225 ~ 36% for 
-D-orbit, and 81/225 ~ 36% for G-orbit. Thus, we can understand the reason why the 5*1, 
Di and Gi orbits in Table ITTl have about 30% occupation probabilities each. 

Figure El^a) demonstrates the radial parts for the 5*1-, Di- and Gi-orbits, the number of 
nodes of which are two, one and zero, respectively. Refiecting the SU(3) character, the radial 
behaviors of the three orbits are similar to those of the harmonic oscillator wave functions 
(unl) with (5 = 4, uo2, U21 and M40, respectively, where (L) denotes the number of nodes 
(orbital angular momentum). We see that the radial parts of the single a-particle orbits 
oscillate widely in the inside region (r < 4 fm). This is due to the strong Pauli blocking 
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effect for the ground state with the compact shell-model-hke structure. The large oscillation 
can also be seen in the reduced width amplitude of the a+^Be{0~^) channel for the ground 
state shown in Fig. El 

Concerning the 0^ state, the occupation probabilities are shown in Table EI We see a 
strong concentration on a single orbit: the occupation probability of the 5*1 orbit is largest, 
amounting to about 70%, and those for other orbits are very small. This means that each 
of the three a particles in the state is in the 5*1 orbit with occupation probability as 
large as about 70%. The radial behavior of the 5*1 orbit is illustrated with the solid line in 
Fig. Efb). We see no nodal behavior but small oscillations in the inner region (r < 4 fm) 
and a long tail up to r ~10 fm. For reference, the radial behavior of the S-wave Gaussian 
function, (fos{r) = A^os(o) exp(— ar^), is drawn with the dashed line in Fig. Iljb), where the 
size parameter a is chosen to be 0.038 fm~^, and NQs{a) denotes the normalization factor. 
The radial behavior of the Si orbit is similar to that of the S'-wave Gaussian function, in 
particular, in the outer region (r > 4 fm), whereas a slight oscillation of the former around 
the latter can be seen in the inner region (r < 4 fm). 

The small oscillation of the Si orbit in the inner region can also be seen in the reduced 
width amplitude of the 0^ state for the a-®Be(0^) channel in Fig. Efa). In order to study 
the origin of the small oscillation, we show in Fig. IHl^b) the results of the reduced width 
amplitudes of the 02" state for the Q;-^Be(0"^) channel, fixing the distance between the 2a 
clusters in ^Be to rQ,a=0.5, 2.5, 4.5 and 6.5 fm. In the case of r^o < 4 fm, we see the nodal 
behavior with the large oscillation in the inner region, coming from the strong Pauli blocking 
effect among the 3a clusters, while the nodal behavior is disappearing and the oscillations 
are getting weaker for the larger r^a (> 4 fm), reflecting the weaker Pauli blocking effect. 
Thus, the small oscillations in the radial behavior of the 5*1 orbit is evidence for the weak 
Pauli blocking effect for the 02" state with the dilute structure. 

The momentum distributions of the a particles, p{k) and k"^ x p{k), are shown for the 
of and 0^ states in Fig. [3 Reflecting the dilute structure of the 0^ state, we see a strong 
concentration of the momentum distribution in the k < 1 fm~^ region, and the behavior of 
p{k) is of the 5-function type, similar to the momentum distribution of the dilute neutral 
atomic condensate states at very low temperature trapped by the external magnetic field . 
On the other hand, the ground state has higher momentum component up to A; ~ 6 fm~^ as 
seen from the behavior of /c^ x p[k) reflecting the compact structure. The above results for 
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the radial behavior of the 5*1 orbit, occupation probabihty and momentum distribution for 
the 0^ state leads us to conclude that this state is similar to an ideal dilute 3a condensate 
with as much as about 70% occupation probability. 

Let us make some remarks on the calculated energy (£'3q,=0.85 MeV) and wave function 
of the 0^ state. They were evaluated under the bound state approximation in the present 
study (see Sec. II). The quite small experimental width for (r=8.5 eV) |3] nieans that the 
bound state approximation is very good to describe the wave function. The complex scaling 



method 



2q is powerful to search for resonant states and calculate the exact energies and 



widths, and is applied easily to the 3a system by slightly modifying the present framework. 
The detailed framework is skipped here and referred to Ref. 2g|. In the present study, we 
investigated the energy of the 0^ state with the complex scaling method. It was found that 
a resonant state, corresponding to the state, appears at -^30=0. 85 MeV with a width less 
than the numerical uncertainty (~100 keV in the present calculation). The results confirm 
that the bound state approximation is good to describe the resonant state. 

It is interesting to compare our results with those given by Matsuura et al. who 
used the normalized spectroscopic amplitude to obtain the bosonic quantities such as the 
single-a orbits and occupation probabilities for the 0^ state in place of the 3a boson wave 
function. According to their results, the occupation probability of the Si orbit {OS orbit 
in Ref. ^^) for the state is about 70%, the value of which is quite similar to ours in 
Table |n] However, the radial behavior of the Si orbit for the 02" state as well as the one 
of the of state given by Matsuura et al. are quite different from our results, and seems 
unnatural. For example, the 5*1 orbit for the 0^ state has as much as 6 ~ 8 nodes and shows 
a behavior similar to that for the Of state, in spite of the fact that the 0^ state has a dilute 
3a condensate structure (see Fig. 6 in Ref. Q|). In addition, the G-orbit for O^" state has a 
prominent peak at r ~ 13 fm, although the state has a shell-model-like compact structure. 
Also the radial behavior of the single-a orbits given in Ref. is hard to understand. This 
may be due to the fact that those authors used the normahzed spectroscopic amplitude in 
place of the 3a boson wave function. 
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B. 2J and 2J states 

The 2t state at E^^^ = -2.83 MeV {E^=4:.U MeV) belongs to the rotational band of 
the ground state starting from the Of state at E^^^ = —7.27 MeV. The calculated energy 
and nuclear radius for 2f in the present study are shown in Table HI E^a = —5.28 MeV 
and 2.45 fm, respectively. The nuclear radius is almost the same as the one for the ground 
state, although the calculated excitation energy is underestimated in comparison with the 



experimental one, in line with what is discussec 



QQH 21, 2^. 



in other papers with the microscopic or 



semi-microscopic 3a cluster model 

The occupation probabilities of the single-a orbits for 2f are demonstrated in Table HTl 
The occupation numbers are concentrated to the first Di and Gi orbits with about 50%. 
Comparing with those for the Of state, we notice the smallness of the occupation number 
for the Si orbit. This feature can be understood from the fact that the 2i state is of the 
SU(3)[/](A/i)j = [444] (04)2+ type with Q=8. The SU(3) state is described as 



I [444] (04)) 2+ = V0.07111|lrf2^), 



+ V0.07111|2slL') - V0A3900\ldlD) - V0.00735|0^1D), 



- V0.00735|lciOG) + ^0.4040810^06"), (56) 

where \{ni){NL)) presents the basis function, \un£{r2a)uNL{ra-2a)) (with 2n + i + 2N + L = 
8), with the harmonic oscillator wave function Une- From the definition of the one-body 
density matrix in Eq. (jHSl), the occupation probabilities for the SU(3) state in Eq. 
are given as 0.07111 for S orbit, 0.07111+0.43900+0.00735=0.51746 for D orbit, and 
0.00735+0.40408=0.41143 for G orbit. Reflecting the character of the SU(3) structure, 
the occupation probability for the orbit in Table |n] is as small as 8.5 %. The radial 
behavior of the single-a orbits, 5*1, Di and Gi ones, is shown in Fig. |Hfa). They are similar 
to those for the 0^ state shown in Fig. IMa)- 

n 

The structure study of ^^C based on the 3a-condensate type wave function |l6[ indicated 
that the 2^ state at E^^^ = 2.6 ± 0.3 MeV with the width of F = 1.0 ± 0.3 MeV Q has a 



structure similar to the O2 state at E^a = 0.38 MeV with the dilute 3a condensation 
The conclusion stems from the result that the 2f state has a large overlap with the single 
condensate wave function of a 3a gas-like structure, the squared value of which amounts 
to about 88%. Thus, it is interesting to study the structure of the 22^ state in the present 
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framework. Since the 23 state is a resonant state with non neghgible width, a continuum 
treatment is requested to estimate exactly the resonant energy and width. 

In order to study the resonant properties of the 2^ state, we take the complex-scahng 
method |26i], which can be apphed easily to the present 3a; system by slightly modifying 
the framework given in Sec. II. The method is powerful to evaluate not only the resonant 
energy and width but also the nuclear radius. The details are again skipped here and we 
refer to Ref. Q. The calculated results are as follows: 1) the 22 resonant state is located at 
i?3Q,=2.3 MeV with r=1.0 MeV, results which are in good agreement with the experimental 
data and 2) the calculated nuclear radius is 4.3 fm, almost the same as that of the 
state. Thus, the 2^ state has a dilute 3a structure. 

It is interesting to study the single-a orbits and occupation probabilities in the 2^ state. 
For this purpose, we need to have the wave function of the 22 state. Since the calculated 
width is not so large in comparison with the resonance energy, the bound state approximation 
is rather good to describe the resonant wave function. The bound state approximation of 
the wave function is obtained within the framework of Sec. II, although the wave function 
gives a large nuclear radius, about 6 fm (see Table H]). Table ITTl illustrates the occupation 
probabilities of the single-a orbits {S-, D- and G-waves) for the 2^ state. We see that the 
occupation probability concentrates on only one orbit, the Di orbit, with occupancy as large 
as 83 %, and the radial behavior of the orbit is likely to be of the D-wave Gaussian-function- 
type with a long tail [see Fig. |H^b)], reflecting a dilute structure. These characteristics are 
quite similar to those for the 0^ state. Thus, we conclude that the 2^ state belongs to the 
3a-condensate structure. 

According to the results in Ref. Q|, it was found that the 2^ state has dominant S'-wave 
between 2a particles and a D-wave between the center-of-mass of the 2a particles and the 
third a. 



This interpretation is consistent with the preset result. The reason is as the follows. From 
the definition of the single-a density matrix in Eq. the single-a density of the 22 state 
is presented as 




) ~ \ue=Q{r2a)UL=2{ra^2a)) ■ 



(57) 



p(r,r') = 3x($(2+)|5(r5'^"^)-r'))(5(ri' 
~ 3 X iV2„ X UL=2{r)UL=2*{r'), 



(cm) 



r) I $(2+)), 



(58) 
(59) 
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where = J dr2aUo* {f2a)uoif2a) ~ 1- Thus, the 2^-state wave function, Eq. (j37jl . has a 
dominant occupation probabihty of the D-orbit, Ul=2- The resuhs are in good agreement 
with the present study. 

C. 3^ state 

The 3~ state at E^^^ =2.37 MeV is an interesting one from the point of view of the 
dilute a condensation. If the state is a condensate with all of the 3a particles in the P 
orbit, there is the possibility of a superfuid with vortex lines, similar to the rotating dilute 
atomic condensate at very low temperature 0. Thus, it is an intriguing problem to study 
the structure in the present framework. 

The calculated energy of the 3^ state is in good agreement with the experimental data 
(see Table H}. The very small width (r*^^'^ = 3.4 keV) indicates that the bound state 
approximation is very good to describe the state. In fact, we checked it theoretically with the 
complex-scaling method, and found that the calculated resonant energy (width) is almost 
the same as the one with the bound state approximation (less than 100 keV, which is the 
numerical uncertainty in the present calculation). Thus, we use the 3~ wave function under 
the bound state approximation to study the characteristics of the state. 

The calculated nuclear radius for the 3~ state is 2.95 fm, the value of which is larger than 
that for the ground state {Of), while it is smaller than that for the 0^ state (see Table H]). 
This suggests that the structure of the 3~ state is intermediate between the shell-model-like 
compact structure (0^") and the dilute 3a; structure (02^). The occupation probability of the 
single-a orbits for the state are shown in Table ITTTl 44.7% for Pi-orbit and 27.9% for Fi- 
orbit. Although the concentration of the single orbit Pi amounts to about 50%, the radial 
behavior of the single-a orbit in Fig. IHl has two nodes in the inner region. However, the 
amplitude of the inner oscillations is significantly smaller than that for the ground state in 
Fig. Efa). The small oscillations indicate the weak Pauli-blo eking effect, and thus, we can 
see the precursor of the 3a condensate state, although the 3~ state is not an ideal rotating 
dilute 3a condensate. 
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D. 1^ state 

The experimental width of the 1^ state at £"3^^=3.57 MeV is as small as r=315 
keV. Q] This means that the bound state approximation is good to describe the state. In 
fact, the calculated energy of the If state under the bound state approximation is -E3a=3.11 
MeV, which is quite similar to that with the complex scaling method (i?3a=3.1 MeV and 
r=0.1 MeV) and in good agreement with the experimental value. 

The calculated nuclear radius, 3.32 fm, is larger than that of the ground state (2.44 fm) 
and the 3j~ state (2.95 fm) but is smaller than that of the one (4.3 fm). The occupation 
probabilities of the a particles in the lj~ state are shown in Table ITTTl 35 % for Pi orbit and 
16 % for F\ orbit. Thus, there is no concentration of the occupation probability to a single 
orbit like the 0^ and 2^ states. Since the a particles in the 1J~ state are distributed over in 
several orbits, the state is not of the dilute a-condensate type. Figure HUl shows the radial 
behavior of the P\ and Fx orbits in the lj~ state. The Pi orbit has two nodes in the inner 
region, the behavior of which is rather similar to the 2P harmonic oscillator wave function. 
However, the Pi orbit has a P-wave Gaussian-type behavior. (Exactly speaking, the orbit 
has one node at the vicinity of the origin, which can not be seen in Fig. E3) Also we see 
the oscillatory behavior of the Pi orbit for < r < 2 fm, similar to the one of the 5*1 orbit 
in the state in Fig. ^b). These interesting behaviors of the Pi orbit indicate some signal 
of the dilute a condensation, reflecting the relatively large nuclear radius (3.32 fm) for the 
\i state. 

E. Structure change of the 0+ state with nuclear radius 

In Sec. IIIA, we found that the state has a dilute 3a structure characterized by the 
nuclear radius as large as about 4.3 fm, in which the a particle occupies in the single orbit 
(5*1) with about 70 % probability, and the radial behavior of the 5*1 orbit is similar to the 
S-wave Gaussian wave function with a very long tail. On the other hand, the Oi" state has 
a compact structure with a nuclear radius of 2.44 fm, where the occupation probabilities 
of the a particles spread out over the S', D and G orbits, amounting to about 30 %, each. 
The feature is much in contrast with that of the 0^ state. The nuclear radius or density 
of ^^C seems to have a close relation with making the compact structure and the dilute "ia 
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structure in the ^^C 0"*" state. Thus, it is very interesting to see the structure change of the 
0"^ state of ^^C by taking the nuclear radius (or density) as parameter. The dependence of 
the occupation probabihties and radial behaviors of the single a-particle orbits in the O"*" 
state on its nuclear radius is investigated with the use of the simple framework given in the 
latter part of Sec. IIC (see Eqs. (jIg|) ~ (IHl|l ). 

Figure ITT] shows the dependence of the energy of ^^C measured from the 3a; threshold on 
the nuclear radius Rn, 2.20 fm < Rn < 4.86 fm, corresponding to a nuclear density 0.15 
^ p/po ^ 1-6 (po denotes the normal density). The energy minimum point appears around 
Rn ~2.4 fm, corresponding to the normal density region. We see the strong repulsion in 
the region of R^ < 2.2 fm, due to the kinetic-energy effect and Pauli-blocking effect, while 
the almost flat region appears at Rn > 4 fm and the energy is positive and small, less than 
1 MeV with respect to the 3a threshold. 

The occupation probabilities of the single-a orbits {Si-, Di-, and Gi-orbits) are shown in 
Fig. ^Jwith respect to the nuclear radius, where denotes the fc-th orbit for the L wave. 
In the region of Rn = 2.2 ~ 2.4 fm (normal density region), the occupation probabilities 
of the a particles spread out over the S, D and G orbits, amounting to about 30 %, each. 
This feature is almost the same as that of the Oj*' state obtained by the 3a OCM calculation, 
the nuclear radius of which is 2.43 fm (see Sec. IIIA). Figure IT^ shows the radial behavior 
of the single-a orbit, 5*1, with respect to the nuclear radius. The ^i-orbit at Rn ~ 2.42 fm 
(Fig. Efa)) has two nodes and the radial behavior is of the 25* harmonic oscillator wave 
function (howf) type, the result of which is almost the same as that of the OJ*" state obtained 
by the 3a OCM calculation (see Fig. Ufa)). Thus, the wave function with Rn ~ 2.4 fm has 
the SU(3)[/](Az/) = [444] (04) character (see Eq. (1^3)). 

Increasing the nuclear radius from Rn = 2.42 fm, the occupation probability of the single- 
a orbits concentrates gradually on a single orbit (^i), and it amounts to be about 90 % at 
i?Ar=4.84 fm (p/po=0.14) in the present calculation (see Fig. IT2|l. The radial behaviors of 
the Si orbit with i?^=2.42, 2.70, 3.11 and 4.84 fm are demonstrated in Figs. CHta), (b), (c) 
and (d), respectively. We can see that increasing the nuclear radius, the internal oscillation 
observed in the orbit with Rn='2A2 fm is gradually disappearing and, finally, the 25- 
type radial behavior transits to the zero-node long-ranged S'-wave type (Gaussian) with the 
occupation probability of about 90 %, approaching an ideal dilute a condensate. The reason 
of why only the 5* wave survives in the case of increasing the nuclear radius is due to the fact 
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that the centrifugal barrier is not at work for the S'-wave a orbit. The S'-wave a particles, 
thus, can move in a nucleus with a given large nuclear radius, although they are confined 

ri 

by the Coulomb potential barrier produced self-consistently [13]. According to the results 
of the 3a OCM calculation (see Sec. IIIA), the a particle in the 0^ state (-Riv=4.3 fm) is 
occupied in the single orbit [S\) with about 70 % probability, the radial behavior of which 
is similar to the S'-wave Gaussian wave function with a very long tail. Their results are 
consistent with those in Figs. ^1 and ^1 
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IV. SUMMARY 



In this work we have investigated the bosonic properties such as single-cc particle orbits 
and occupation numbers in the J'^=0+, 2+, 1", and 3~ states of ^^C around the 3a thresh- 
old within the framework of the 3a OCM (orthogonality condition model). The 3a OCM 
equation is based on the equation of motion for the Na bosons derived from the microscopic 
Na cluster model theory. The experimental energy spectra for Of, Of , 2f, 1^, and 3^^ are 
reproduced well with the So; OCM. 

The main results to be emphasized here are as follows. 

(1) The O2" state at £^3^^=0.38 MeV has a dilute 3a structure characterized by the nu- 
clear radius as large as about 4.3 fm. The analysis of the single-a orbits and occupation 
probabilities for the dilute state shows that the a particle is occupied in a single orbit (Si) 
with about 70 % probability, and the radial behavior of the Si orbit is similar to the S-wave 
Gaussian wave function with a very long tail. The momentum distribution of the a particle 
illustrates the (5-function like behavior, similar to the momentum distribution of a dilute 
neutral atomic condensate states at very low temperature, a feature which eventually can 
be measured experimentally. These results give more theoretical evidence that the O2 state 
is a dilute 3a condensate. On the other hand, the Of state has a compact structure with a 
nuclear radius of 2.44 fm. The occupation probabilities of the a particles spread out over the 
S, D and G orbits, amounting to about 30 %, each, the results of which comes from the fact 
that the Oi" state is characterized by the nuclear SU(3) wave function, [/](A//) = [444] (04). 
The feature is much in contrast with that of the Of state. 

(2) In order to understand further the characteristic structure of the two 0+ states, we have 
studied the single-a orbital behavior in the ^'^C{0~^) state by taking the nuclear radius 

(or density p/po) as parameter, 2.42 < Rj^ < 4.84 fm, (0.15 < p/ Po < 1.2), where po denotes 
the normal density). We found that the single-a orbits in the ^^C(0+) state are smoothly 
changed with the nuclear radius Rn, and their behavior is classified into the following three 
types, depending on the value of it! at: i) at Rn ~ 2.4 fm (p/po ~ !)> we have two-nodal S- 
orbit (25"), one-nodal D-orbit (ID) and zero-nodal G-orbit (OG) with about 30% occupation 
probability, each, characterized by a nuclear SU(3) wave function, ii) increasing the nuclear 
radius from i?jv ~ 2.4 fm, the occupation probability of the single-a orbits concentrates 
gradually on a single ^-orbit, in which the two-nodal behavior is disapperaing, and then. 
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iii) at Rn ~ 4 fm (p/po ~ 0.2), there appears a dominant zero-nodal Gaussian (05'-type) 
orbit with a very long tail, the radial behavior of which is similar to that of the O2 state 
in ^^C as mentioned above. The structure change is caused mainly by the Pauli-blocking 

effect, the strength of which depends dominantly on the nuclear radius in the present 
framework. 

(3) The structure of the 2^ state at £^|^^=2.6 ± 0.3 MeV with T = 1.0 ± 0.3 MeV 
was studied with the present 3a OCM and the complex-scahng method. We found that 
the 2^ resonant state appears at £^3^=2.3 MeV with F =1.0 MeV, in agreement with the 
experimental data, and the calculated nuclear radius is 4.3 fm, similar to that of the O2 state. 
The 2^ wave function obtained with the 3a OCM was used to study the bosonic properties 
of the state. It was found that the occupation probability of the a particle concentrates 
only on the Di orbit, amounting to be as large as about 80 %, and the radial behavior is 
of the D-wme Gaussian type with long tail. The characteristics of the boson properties 
in 22 is quite similar to those in O2 at £'3^^=0.38 MeV. Thus, the 22 state has the dilute 
3Q;-condensate-like structure. On the other hand, the 2^ state has a compact structure 
with the nuclear radius, 2.44 fm, like the ground state. The occupation probabilities of 
the a particles spread out over the D and G orbits, amounting to about 56 % and 33 %, 
respectively, reflecting the SU(3) character of the 2i state. 

(4) We investigated the a bosonic properties of the negative parity states, 1 j" at ^3^^=3.57 
MeV and 3^ at £;|^*'=2.37 MeV. Their nuclear radii are 3.32 and 2.95 fm, respectively, which 
are larger than that of the ground state (O^") but smaller than that of O2". The calculated 
occupation probabihties of the a particles in those states show that there is no concentration 
on a single a orbit like in the 0^ and 2^ states. The results indicates that the 1^ and 3J~ 
states are not of the dilute 3a condensate. The radial behavior of the P- and F-wave single 
a orbits, however, suggests that small components of the 3a condensation exist even in the 
negative parity states, which is reflected by their relatively large nuclear radii. 
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TABLE I: Calculated energies (E^a) and nuclear radii {y {r%)) for the 0+, 2+, 3~ and 1~ states 
of ^^C together with the a-a and a-2a rms distances (V {faa) and ^ {i^a-2a))- "^^^ energy E^ce is 
one measured from the 3a threshold. The values in parentheses denote the experimental ones. All 
energies and nuclear radii (rms distances) are given in units of MeV and fm, respectively. 





jpcal 










or 


-7.27 


(-7.27) 


2.44 


3.02 


2.61 


0+ 


0.85 


( 0.38) 


4.31 


6.84 


5.93 


2^ 


-5.28 


(-2.83) 


2.45 


2.94 


2.55 


2^ 


2.43* 


( 2.6) 


6.12* 


10.2 


8.80 


3r 


1.52 


( 2.37) 


2.96 


4.10 


3.56 


ir 


3.11 


( 3.57) 


3.32 


4.87 


4.23 



*According to the complex-scaling method, the resonant energy and width of the 2^ state are E^a = 2.3 
MeV and T = 1.0 MeV, respectively, with ^(r^)=4.3 fm. See text. 
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TABLE II: Occupation numbers of the k-th a-orbits with S-, D- and G-waves for the 0^ and 2+ 
states of ^^C obtained by diagonalizing the one-body density matrix in Eq. (|33() . The values in 
parentheses denote the occupation probabihties. 





k 


S 






G 




1 


1.05 (35.0 


%) 


1.06 (35.3 %) 


0.82 (27.3 %) 




Z 


U.Ui U.o 


/o J 


U.Ui 1^ u.u /o) 


n nn (' n n 
U.UU 1^ u.u /o j 




3 


0.00 ( 0.0 


%) 


0.01 ( 0.0 %) 


0.00 ( 0.0 %) 




total 


1.07 (35.6 


%) 


1.07 (35.6 %) 


0.82 (27.3 %) 




1 


2.16 (72.0 


%) 


0.19 ( 6.3 %) 


0.08 ( 2.7 %) 




9 


U.ZU 0. / 


/o ) 


U.UD z.u /O J 


n nfi i' 9 n cw^i 
U.UD z.u /o) 




3 


0.02 ( 0.7 


%) 


0.02 ( 0.7 %) 


0.01 ( 0.3 %) 




total 


2.38 (79.3 


%) 


0.29 ( 1.0 %) 


0.16 ( 5.3 %) 




1 


0.25 ( 8.5 


%) 


1.69 (56.2 %) 


1.00 (33.3 %) 




Z 


U.UU u.u 


/o ) 


n no ( n 7 

U.KJZ [U.I /o) 


n nn ( rt n 
U.UU u.u /o J 




3 


0.00 ( 0.0 


%) 


0.00 ( 0.0 %) 


0.00 ( 0.0 %) 




total 


0.26 ( 8.7 


%) 


1.71 (57.0 %) 


1.00 (33.3 %) 




1 


0.31 (10.3 


%) 


2.50 (83.3 %) 


0.05 ( 1.7 %) 




2 


0.02 ( 0.7 


%) 


0.03 ( 1.0 %) 


0.00 ( 0.0 %) 




3 


0.00 ( 0.0 


%) 


0.01 ( 0.3 %) 


0.00 ( 0.0 %) 




total 


0.33 (11.0 


%) 


2.56 (85.3 %) 


0.06 ( 2.0 %) 
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TABLE III: Occupation numbers of the k-th a-orbits with P- and F-waves for the 3^ and 1~ 
states of ^^C obtained by diagonalizing the one-body density matrix in Eq. (|33() . The values in 
parentheses denote the occupation probabihties. 





k 


P 






1 


1.34 (44.7 %) 


0.84 (27.9 %) 







U. iZ 4.U /O j 


U.Zo 1^ (.0 /o j 




3 


0.06 ( 1.9 %) 


0.02 ( 0.8 %) 




total 


1.54 (51.4 %) 


1.09 (36.4 %) 




1 


1.06 (35.3 %) 


0.47 (15.8 %) 




2 


0.53 (17.8 %) 


0.26 ( 8.6 %) 




3 


0.08 ( 2.6 %) 


0.08 ( 2.6 %) 




total 


1.75 (58.5 %) 


0.84 (28.1 %) 
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(a) (b) (c) 



FIG. 1: Jacobian coordinate systems for the 3a model of C. The three Jacobian coordinates, 
(a), (b) and (c), correspond respectively to the 3a relative wave functions, <I>j"(12,3), <^j°(23, 1) 
and ^f'{3l,2) in Eq. 
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FIG. 2: Density distribution of the a particle for the Of (sohd hne) and (dotted) states. 
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FIG. 3: Correlation functions, (a) faa and (b) fa-2a, for the 0^ (solid line) and 0^ (dotted) 
states. 
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FIG. 4: Radial parts of the single a orbits, (a) Si (solid line), Di (dashed) and Gi (dotted), in 
the Oi state, and (b) the Si (solid) orbit in the state and S'-wave Gaussian function (dotted), 
rtpos, with the size parameter a = 0.038 fm~^ (see text). Note that all of the radial parts in figures 
are multiplied by r. 36 
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FIG. 7: Momentum distribution of the a particle, (a) and (b) k'^p{k), for the 0^ (sohd Hne) 
and O2 (dotted) states. 



39 
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FIG. 8: Radial parts of the single a orbits, (a) 5*1 (solid line), D\ (dashed) and G\ (dotted), in the 
2^ state, and (b) the D\ (solid) orbit in the 22 state. Note that all of the radial parts in figures 
are multiplied by r. 
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FIG. 9: Radial parts of the single a orbits, Pi (solid line) and F\ (dotted), in the 3^ state. Note 
that all of the radial parts in figures are multiplied by r. 
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FIG. 10: Radial parts of the single a orbits, Pi (solid line) and Fi (dotted), in the 1^^ state. Note 
that all of the radial parts in figures are multiplied by r. 
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FIG. 11: Dependence of the energy of the ^^0(0"'") state measured from the 3a threshold on its 
nuclear radius. 
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FIG. 12: Dependence of the occupation probabilities of the single-a orbits (-Si-, Di- and L'l-orbits) 
in the ^^0(0"*") state on its nuclear radius. The sohd (dotted and dot-dashed) hne corresponds to 
the S'l-orbit {Di- and Gi orbits, respectively). 
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FIG. 13: Radial behaviors of the Si-orbit in the ^^0(0+) state with (a) Rn=2A2 fm, (b) i?Ar=2.70 
fm, (c) Rn=3.11 fm, and (d) i?Ar=4.84 fm, where Rn denotes the nuclear radius of the ^^C(0"*") 
state. 
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